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Superconducting BCS versus Fulde-Ferrell-Larkin-Ovchinnikov states of quasiparticles 
with spin dependent mass and their distinguishability 
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The first observation of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconducting state and 
a subsequent detection of the spin-dependent effective masses of quasiparticles in the CeCoIns heavy 
fermion system are combined into a single theoretical framework. The appearance of the spin-split 
masses extends essentially the regime of temperatures and applied magnetic fields, in which FFLO is 
observable and thus is claimed to be very important for the FFLO detectability. We also stress that 
the quasiparticles composing Cooper pair become distinguishable in the nonzero field. The analysis 
is performed within the Kondo-lattice limit of the finite-C/ Anderson - lattice model containing both 
the mass renormalization and real-space pairing within a single scheme. 

PACS numbers: 74.70.Tx, 71.27. +a, 74.20.Mn 
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The superconductors discovered at the beginning of 
this decade are termed as unconventional because the 
breakdown of the basic symmetries such as the spatial [lj 
or time [2[ or the spin singlet-triplet Cooper pair mixing 
3 1 have been discovered in the heavy-fermion compounds 
J|. These features alone lead to a nonstandard behav- 
ior even when the Bogolyubov-de Gennes quasiparticle 
approach is applicable. Such an approach is based on 
the concept of Landau-Fermi liquid, albeit almost local- 
ized [5[. A separate, largely unanswered question, is con- 
cerned with the pairing in the non-Fermi (non-Landau) 
liquids ((J and particularly, in high temperature super- 
conductors 0]. In this situation, it seems proper to se- 
lect cases, for which a limited number of novel factors 
such as the spin-dependence of quasiparticle mass [8[ and 
the real-space pairing induced by strong correlations @, 
can be incorporated into an effective Fermi-liquid picture, 
and thus become tractable theoretically [10| . We suggest, 
that the heavy-fermion systems such as CeCoIn 5 [ll[ and 
the specific organic compounds [11] can be selected as 
such test cases. 

In this paper we take into account spin dependent 
quasiparticle masses @ observed recently [l3[ and ana- 
lyze in detail the relative stability of the Bardeen-Cooper- 
Schrieffer (BCS) phase against the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state [l5| in an applied mag- 
netic field H and at temperature T > 0. The basic 
question we tackle, apart from a detailed quantitative 
analysis, is why the FFLO phase became so easily ob- 
servable in this strongly correlated system, after almost 
forty years of experimental search for that phase exis- 
tence. By introducing the spin-split masses m a (a = ±1 
is the quasiparticle spin quantum number), one deals, 
in our view, with a fundamentally new situation, as by 
switching on the applied field one may transform the sys- 
tem of quantum mechanically indistinguishable quasipar- 



ticles into their distinguishable correspondants. This cir- 
cumstance produces not only changes of the normal-state 
properties such as the appearance of itinerant-electron 
metamagnetism [l6| , but also leads to the essential mod- 
ification of the single Cooper-pair bound state composed 
of two distinguishable particles [17f . In this manner, the 
spin-dependent renormalization of the quasiparticle mass 
introduces by itself a non-Landau feature to fermionic liq- 
uid and appears only when the correlated electrons are 
treated within a non-perturbational approach 0, [lfjl . [f8| . 
Here we concentrate first on the condensed state of the 
pairs, which retain their indistinguishability and discuss 
subsequently the properties of a single Cooper pair, with 
partners becoming distinguishable in the quantum me- 
chanical sense, in the applied field. 

We start with the effective Anderson-lattice model in 
the large but finite- U limit [HI], which has the form 
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The first term represents the conduction (c) electrons, 
the second the originally localized (/) electrons with the 
site double occupancies projected out {Nia = f^fia), the 
third - the projected f - c hybridization, and the fourth 
the real-space pairing part in the leading order, with the 
projected spin-singlet pairing operators 



&L = ^ [^tC 1 - N n)cml - a.V 1 - %)cmr]- (2) 

Note that the pairing part disappears if the /-/ in- 
traatomic interaction U — > oo. The crucial features of 
([T]) are the inclusion of the projection which leads to the 
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spin-dependent renormalization of the hybridization and 
in turn, the spin-split quasiparticle masses, as well as 
of the virtual (f-c) hopping induced real-space pairing, 
both within a single scheme. In the Kondo-lattice limit 
i.e. when n f = (£ CT iVi CT (l " N ia )) = 1-5, with S < 1) 
and in the saddle-point or Gutzwiller approximations || , 
the Hamiltonian can be brought (up to a constant 
term) into a single-narrow band form with the BCS-type 
of pairing Q including the Zeeman term (field H ^ 0) 
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where the operators ip^ ~ /jL, q a ^u being the band 
energy of itinerant / electrons induced by the hybridiza- 
tion, 7k is a complicated function [9( of quasimomen- 
tum k, and Vkk' = -4V£V£/(ej(ef + u )) is tne strength 
of the pairing potential, taken in the following for in- 
traatomic form of hybridization, Vk = V. The factor 
q a ~ (1 — n/)/(l — nf a ) is the hybridization renormal- 
ization implying the spin dependent mass renormaliza- 
tion, m a /mo = (Q x . Physically, the pairing results 
from the Kondo-type self-screening of the heavy quasi- 
particles. The renormalization factor of the pairing part 
Ra-a ~ (laqw) 1 ^ 2 will be regarded as a constant reduc- 
ing the pairing potential magnitude V$ = ^ (^+u) ' ^ n 
what follows we take the form of 7k as the d— wave in 
the two-dimensional limit. 

A few important physical remarks conveying principal 
features of our approach are in place here. First, the 
original atomic / electrons acquire band properties by a 
three-step process: hopping / — > c from the atomic to the 
bare conduction band states, followed by a propagation 
in the conduction band, and a subsequent deexcitation 
c — > /. Hence, the hopping amplitude for / states is of 
the order (V/ef) 2 t mn . Second, the mechanism of pairing 
is of analogous origin as that in t-J model with a different 
renormalization factors for the pairing part [19j |. Third, 
the strongly correlated and hybridized electrons form an 
almost localized Fermi liquid fioj of electrons in a very 
narrow band of / states. The last assumption means that 
the number of quasiparticles characterized by rif (related 
to e. g. cerium valency by Ce +4_n/ ) is constant as a func- 
tion of T and H which is valid only at low temperature 
T <T K [T K is the effective Kondo temperature [||). Fi- 
nally, Tk plays the role of the /-electron bandwidth and 
its value is of the order 2z\t\(V/e f ) 2 (l-n f ) ~ lCT 2 t = zl, 
where t is the amplitude of the hopping t mn between z 
nearest neighbors. Since t is a fraction of eV, then Tk 
is in the regime 10 — 100 K. This is the reason why the 
heavy electrons are so sensitive to Zeeman fields of the 
order of 10 T. 

We discuss now the relative stability of the BCS and 
the FFLO condensed states, the latter defined by the 
center-of-mass momentum Q ^ of the Cooper pairs and 



by the gap amplitude Aq . An essential factor is that the 
states for H ^ have spin-dependent masses m a , as the 
field acts on both quasiparticles resulting from the broken 
pairs, as well as from gapless regions of Ak = Ak,Q. 
Having in mind different gap symmetries, we assume that 
Ak = Acj7(k), where the symmetry expressed via 7(k) 
is superposed on Q dependence of the gap magnitude. 
We analyze explicitly the case of quasi-two-dimensional 
superconductor [20| in a clean limit, where Pauli-limiting 
situation [2l| is well defined. Also, as we consider real 
space pairing caused by strong correlations, we require 
a detailed dispersion relation in the full band, which is 
parametrized as follows 

£k(T = q<y[— 2t(cosk x + cos k y ) + At 1 cos k x cos k y ] , (4) 

with t and t' being hopping integrals between the 
nearest and the next-nearest neighbors (we take t' jt = 
0.5). This electronic structure represents properly the 
situation in both CeCoIns [22j and of k — (BEDT — 
TTF)2Cu(NCS)2- Using the standard diagonalization 
techniques we obtain the quasiparticle energies in the su- 
perconducting state in the form similar to that in [20l[23|. 



£k,Q,± - + ^(^T ~ e -k+Ql) 

± \{{e^ + 6-k+Qi - 2 M ) 2 + 4|A k , Q | 2 ]^ (5) 

where fj, is the chemical potential and the gap is deter- 
mined from the self-consistent equation 
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This equation must be supplemented by the correspond- 
ing equations for /i and n a . In effect, one has to minimize 
the free-energy functional T — T(T, H; Aq, Q, n a ) of the 
form 



T = ~k B Tj2 Ml + exp(-£;k,Q, a /A: B T)] + fiNn f 

k a=± 

+ ]>>- k -Qi " ^ " 9»bH) + N\A Q \ 2 /V , (7) 



(N being the number of sites), with respect to Q, 
Aq, and the band fillings n a (subject to the constraint 
n a + n-g- = nf). The solutions were determined by the 
two methods to check for the accuracy: first by sampling 
the Brillouin zone and second by integrating over bare 
density of states /?o(e) with a variable grid depending on 
the steepness of po(e) behavior, particularly near the van 
Hove singularity. 

In Fig. lab we display the field dependence of the mass 
enhancement in the normal state close to the integer nf. 
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FIG. 1: (Color online). Spin-dependent mass enhancement 
vs. H for selected values of both the filling n/ and tempera- 
ture T. The effect is particularly strong for n/ very close to 
unity and for T < T K - 



Note a strong temperature dependence associated with 
the reduction of the magnetization (n-f — fij.). Taking t = 
20 K, we estimate the regime of physical fields H < 20 T 
as limited by figH/t < 0.1, where the mass splitting is 
indeed essential. 

To characterize the stability of the BCS and FFLO 
states, we plot the phase diagrams in Fig. 2ab with spin- 
independent (average m av = {fn-\ + mjj/2) and with the 
spin-dependent {m a ^ m^) masses, respectively. One 
notes the essential extension of the stability regime (of 
the FFLO phase) on the H — T plane in the latter case. 
This is the reason, we claim, why the FFLO state has 
been observed for the first time in CeCoIns, where the 
spin-dependent masses were observed fl3T ] . For complete- 
ness, in Fig. 3ab we plot the values of the gap Aq and the 
optimal Q values, respectively (Aq =0 in the BCS state 
and Ao = A then). The order parameter A jumps across 
the BCS-FFLO border signalling the lst-order transition. 
The transition to the normal phase seems to be quasicon- 
tinuous. The dashed lines in Fig. 2ab mark the critical 
field H C 2 for the BCS state. Also, the uppermost critical 
field is spectacularly curved upward only if we include 
the spin mass splitting, which is observed in the organic 
system [231 ], as well as the discontinuous nature of the 
BCS-FFLO transition. 

One should note that the free-energy differences at the 
transition are of the order of 10~ 4 t at most, so the ac- 
curacy of extracting the solution of the self-consistent 
integral equations for Aq, n CT , and fi is particularly in- 
tricate. This circumstance is also the reason behind the 
statement that within our numerical accuracy, we cannot 
decisively state that the FFLO-normal state phase tran- 
sition is continuous. The last remark is illustrated in Fig. 
3ab, where we plot the gap amplitude on the H — T plane. 
In Fig. 3b we provide the value of Q, which is a well de- 
fined quantity in FFLO regime (Q x = Q y w 0.2 — 0.8 in 
7r/a units, where a is the lattice constant). 
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FIG. 2: (Color online). Phase boundaries for a d-wave su- 
perconductor with both the spin-independent masses (a) and 
with the spin-split masses (b). The FFLO-BCS transition line 
is discontinuous. The dashed line marks the stability limit of 
the BCS state as determined by the value of the critical field 
H C 2- The values of parameters: n/ = 0.97 and Vo = 10 K. 



From the results above it follows that the spin- 
dependent masses play an important quantitative role 
in the stability of the phases. They lead also to a qual- 
itative, if not fundamental, change of the description of 
individual particles (cf. Fig. 1) and of the single-Cooper 
pair state. Namely, since the mass splitting ro-f — mi ~ 
( n /T — n /l)i the indistinguishable particles in the field 
(or spin-polarization) absence, become distinguishable for 
H > 0. To demonstrate this we note first, that since the 
masses are labeled by the spin quantum number, this 
means that the spin part of the wave function should be 
chosen either as xt(1)xi(2) or as xj(l)xf (2), not in the 
usual singlet form [x T (l)x;(2) - X; (1) Xt (2)]/V2. This is 
what we mean [251 ] by the quasiparticle distinguishability. 
To test our concepts we reformulated 0, [2y] the orig- 
inal Cooper problem including the momentum Q =/= 0. 
The basic Cooper-pair properties are displayed in Fig. 
4ab. In Fig. 4a we plot the pair binding energy A in 
a few situations: the uppermost solid curve is the de- 
pendence Aq(H) with the spin dependent masses, start- 
ing from m av — 304 m and Q ^ 0; the solid line just 
below reflects the corresponding dependence for Q = 0. 
The modified spin- wave function essentially increases the 
binding energy [l6j]. In Fig. 4b we plot the value of |Q| 
vs. H. The jump to the value Q ~ \kF] — ^fx! (fcpv 
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FIG. 3: (Color online). Magnitude A of the order parameter 
on the H — T plane (a) and the center-of-mass momentum 
Q (b) for the same situation as above in the case with spin- 
dependent masses. 



denotes the Fermi wavevector for the c-subband) is rep- 
resented by the empty circle in Fig. 4a. The solid circles 
define the field values destabilizing the pair bound state 
through the spin flip. The stable state of the single-pair 
state with Q ^ in H > 2.5Tis thus the prerequisite of 
the FFLO state. 

In summary, we have applied the concept of the spin- 
direction dependent mass for a fermionic liquid composed 
of quasiparticles that condense into the BCS or FFLO 
states. This basic property extends the regime of the 
FFLO phase stability on the expense of the BCS state. 
The FFLO regime is stabilized by the circumstance that 
the mass difference enhances behavior of system corre- 
sponding difference in the density of states p a (e). We 
also point out to the basic nontrivial feature of the quasi- 
particles composing Cooper pair. Namely, we can study 
an evolution of the two-particle state from the limit of 
indistinguishable to distinguishable quantum particles as 



a function of physically controllable parameter. Such re- 
sult appears in both the Gutzwiller (or mean field) and 
the dynamic mean-field levels [1, [13] • The situation close 
to the metamagnetic point requires a separate analysis of 
an effective field appearing on the saddle-point level (cf. 

1). 
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FIG. 4: (Color online), (a) field dependence of single-pair 
binding energy with spin-split masses (Am 7^ 0) and without 
(Am = 0), for both center of mass momentum Q 7^ and 
Q = 0. The empty circle marks the transition to the state 
with Q roughly equal to the difference of the Fermi momenta 
fci?T an d fl»Fii whereas those on the abcissa are the Pauli fields 
for the bound-state destruction, (b) Q vs. H for the stable 
pair state (the spin wave-function symmetry is marked). 
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